Observation of the decay ip(2S) K* (892) K + c.c. 

M. Ablikim 1 , J. Z. Bai 1 , Y. Ban 10 , J. G. Bian 1 , X. Cai 1 , J. F. Chang 1 , H. F. Chen 16 , 
H. S. Chen 1 , H. X. Chen 1 , J. C. Chen 1 , Jin Chen 1 , Jun Chen 6 , M. L. Chen 1 , Y. B. Chen 1 , 

S. P. Chi 2 , Y. P. Chu 1 , X. Z. Cui 1 , H. L. Dai 1 , Y. S. Dai 18 , Z. Y. Deng 1 , L. Y. Dong 1 , 
S. X. Du 1 , Z. Z. Du 1 , J. Fang 1 , S. S. Fang 2 , C. D. Fu 1 , H. Y. Fu 1 , C. S. Gao 1 , Y. N. Gao 14 t, 
M. Y. Gong 1 , W. X. Gong 1 , S. D. Gu 1 , Y. N. Guo 1 , Y. Q. Guo 1 , Z. J. Guo 15 , F. A. Harris 15 , 
K. L. He 1 , M. He 11 , X. He 1 , Y. K. Heng 1 , H. M. Hu 1 , T. Hu 1 , G. S. Huang 1 * , L. Huang 6 , 
X. P. Huang 1 , X. B. Ji 1 , Q. Y. Jia 10 , C. H. Jiang 1 , X. S. Jiang 1 , D. P. Jin 1 , S. Jin 1 , Y. Jin 1 , 
Y. F. Lai 1 , F. Li 1 , G. Li 1 , H. H. Li 1 , J. Li 1 , J. C. Li 1 , Q. J. Li 1 , R. B. Li 1 , R. Y. Li 1 , S. M. Li 1 , 
(N| ' W. G. Li 1 , X. L. Li 7 , X. Q. Li 9 , X. S. Li 14 , Y. F. Liang 13 , H. B. Liao 5 , C. X. Liu 1 , F. Liu 5 , 

Fang Liu 16 , H. M. Liu 1 , J. B. Liu 1 , J. P. Liu 17 , R. G. Liu 1 , Z. A. Liu 1 , Z. X. Liu 1 , F. Lu\ 
G. R. Lu 4 , J. G. Lu 1 , C. L. Luo 8 , X. L. Luo 1 , F. C. Ma 7 , J. M. Ma 1 , L. L. Ma 11 , Q. M. Ma 1 , 

X. Y. Ma 1 , Z. P. Mao 1 , X. H. Mo 1 , J. Nie 1 , Z. D. Nie 1 , S. L. Olsen 15 , H. P. Peng 16 , 
N. D. Qi 1 , C. D. Qian 12 , H. Qin 8 , J. F. Qiu 1 , Z. Y. Ren 1 , G. Rong 1 , L. Y. Shan 1 , L. Shang 1 , 
D. L. Shen 1 , X. Y. Shen 1 , H. Y. Sheng 1 , F. Shi 1 , X. Shi 10 , H. S. Sun 1 , S. S. Sun 16 , Y. Z. Sun 1 , 
Z. J. Sun 1 , X. Tang 1 , N. Tao 16 , Y. R. Tian 14 , G. L. Tong 1 , G. S. Varner 15 , D. Y. Wang 1 , 
* ; J. Z. Wang 1 , K. Wang 16 , L. Wang 1 , L. S. Wang 1 , M. Wang 1 , P. Wang 1 , P. L. Wang 1 , 

p; S. Z. Wang 1 , W. F. Wang 1 , Y. F. Wang 1 , Zhe Wang 1 , Z. Wang 1 , Zheng Wang 1 , 

Z. Y. Wang 1 , C. L. Wei 1 , D. H. Wei 3 , N. Wu 1 , Y. M. Wu 1 , X. M. Xia 1 , X. X. Xie 1 , B. Xin 7 , 
G. F. Xu 1 , H. Xu 1 , Y. Xu 1 , S. T. Xue 1 , M. L. Yan 16 , F. Yang 9 , H. X. Yang 1 , J. Yang 16 , 
S. D. Yang 1 , Y. X. Yang 3 , M. Ye 1 , M. H. Ye 2 , Y. X. Ye 16 , L. H. Yi 6 , Z. Y. Yi 1 , C. S. Yu 1 , 
G. W. Yu 1 , C. Z. Yuan 1 , J. M. Yuan 1 , Y. Yuan 1 , Q. Yue 1 , S. L. Zang 1 , Yu Zeng\Y. Zeng 6 , 
B. X. Zhang 1 , B. Y. Zhang 1 , C. C. Zhang 1 , D. H. Zhang 1 , H. Y. Zhang 1 , J. Zhang 1 , 
J. Y. Zhang 1 , J. W. Zhang 1 , L. S. Zhang 1 , Q. J. Zhang 1 , S. Q. Zhang 1 , X. M. Zhang 1 , 
X. Y. Zhang 11 , Y. J. Zhang 10 , Y. Y. Zhang 1 , Yiyun Zhang 13 , Z. P. Zhang 16 , Z. Q. Zhang 4 , 
D. X. Zhao 1 , J. B. Zhao 1 , J. W. Zhao 1 , M. G. Zhao 9 , P. P. Zhao 1 , W. R. Zhao 1 , X. J. Zhao 1 , 

Y. B. Zhao 1 , Z. G. Zhao 1 *, H. Q. Zheng 10 , J. P. Zheng 1 , L. S. Zheng 1 , Z. P. Zheng 1 , 
X. C. Zhong 1 , B. Q. Zhou 1 , G. M. Zhou 1 , L. Zhou 1 , N. F. Zhou 1 , K. J. Zhu 1 , Q. M. Zhu 1 , 
Y. C. Zhu 1 , Y. S. Zhu 1 , Yingchun Zhu 1 , Z. A. Zhu 1 , B. A. Zhuang 1 , B. S. Zou 1 . 

(BES Collaboration) 

1 Institute of High Energy Physics, Beijing 100039, People's Republic of China 
2 China Center for Advanced Science and Technology (CCAST), Beijing 100080, People's 

Republic of China 

3 Guangxi Normal University, Guilin 541004, People's Republic of China 



Preprint submitted to Elsevier Science 



7 February 2008 



16 



4 Henan Normal University, Xinxiang 453002, People's Republic of China 
5 Huazhong Normal University, Wuhan 430079, People's Republic of China 
6 Hunan University, Changsha 410082, People's Republic of China 
7 Liaoning University, Shenyang 110036, People's Republic of China 
8 Nanjing Normal University, Nanjing 210097, People's Republic of China 

9 Nankai University, Tianjin 300071, People's Republic of China 

10 Peking University, Beijing 100871, People's Republic of China 

11 Shandong University, Jinan 250100, People's Republic of China 

12 Shanghai Jiaotong University, Shanghai 200030, People's Republic of China 

13 Sichuan University, Chengdu 610064, People's Republic of China 

14 Tsinghua University, Beijing 100084, People's Republic of China 

15 University of Hawaii, Honolulu, Hawaii 96822, USA 
University of Science and Technology of China, Hefei 230026, People's Republic of China 
17 Wuhan University, Wuhan 430072, People's Republic of China 
18 Zhejiang University, Hangzhou 310028, People's Republic of China 



* Current address: University of Michigan, Ann Arbor, Michigan 48109, USA 
"I" Current address: Purdue University, West Lafayette, Indiana 47907, USA. 



Abstract 



Using 14 million ip(2S) events collected with the BESII detector, branching fractions of ip(2S) — > K* (892) + 7\ +c.c. 
and K*(892)°K° + c.c. are determined to be: B(ip(2S) -> K* (892) + K~ + c.c.) = (2.91^ ± 0.4) x 1(T 5 and 
B(ip(2S) -> K* (892) c 'K° + c.c.) = (13.3±^ ± 1.7) x 10" 5 . The results confirm the violation of the "12% rule" for 
these two decay channels. A large isospin violation between the charged and neutral modes is observed. 
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1. Introduction 

One of the longstanding mysteries in heavy quarko- 
nium physics is the strong suppression of -0(25') decays 
to the vector-pseudoscalar (VP) meson final states, pn 
and K*(892)K + c.c, referred to as the "p7r puzzle". 
In perturbative QCD (pQCD), hadronic decays of the 
J ftp and ip(2S) are expected to proceed dominantly 
via three gluons or a single direct photon, with widths 
proportional to the square of the cc wave function at 
the origin, which is well determined from dilepton de- 
cays. Thus, it is reasonable to expect, for any hadronic 
final state h, the J/ip and ip(2S) decay branching frac- 
tions should satisfy the so-called "12% rule"(|l|) 



Qh 



gMgg) -> h) _ B(iJj(2S) -> e+e") 
B(J/ip -> h) ~ B{J/ijj -> e+e~) 



12%, 



where the leptonic branching fractions are taken from 
the Particle Data Group (PDG) tables (0). It was 



first observed by the Markll experiment that, while 
this rule works reasonably well for a number of exclu- 
sive hadronic decay channels, it is severely violated 
for the vector-pseudoscalar meson (VP) final states, 



pn and K*(892) + K~ + c.c. p. Preliminary BESI re- 
sults confirm the Markll measurements at higher sen- 
sitivity (4). This anomaly has generated much interest 
and led to a number of theoretical explanations 
More precise experimental results are required to dis- 
tinguish between them. 

In this paper, the branching fractions of charged 
and neutral tp(2S) — * K*(892)K + c.c. decays are 
reported, based on a sample of 14.0 x 10 6 (1 ± 4%) 
■0(25) events Jjj) collected with the Beijing Spectrom- 
eter (BESII) (|3) at the Beijing Electron-Positron Col- 
lider (BEPC). 

In this analysis, a GEANT3 based Monte Carlo 
package (SIMBES) with detailed consideration of the 
detector performance (such as dead electronic chan- 
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nels) is used. The consistency between data and Monte 
Carlo has been carefully checked in many high pu- 
rity physics channels, and the agreement is reason- 
able Q). The generator KSTARK Q, which simu- 
lates ip(2S) -> K*(892)K + c.c. events, together with 
SIMBES, is used to determine detection efficiencies. 

2. Event selection 

Candidate events for this decay mode have the final 
state K^K ± — > 7r + TT K^nT . They are required to 
satisfy the following general selection criteria: (i) The 
number of charged particles must be equal to four with 
net charge zero, (ii) Each charged track is required to 
be well fitted to a three dimensional helix and be in 
the polar angle region | cos#| < 0.8. (iii) Background 
from ip(2S) — > 7r + 7r~ J/ip, J/ip — > X is removed by 
requiring that the recoil mass of any two oppositely 
charged tracks satisfies 

KZcott = V(Ecra -E+- - (p + +p_)2 

£ (3.05,3.15) GeV/c 2 , 

where E + (E_) and are the assumed tt + (tt^) 

energy and momentum, respectively. 

The Kg is identified through the decay Kg — > 
7r + 7r~. The intersections of all pairs of oppositely 
charged tracks, assumed to be pions, are found as 
described in Ref . jlll) . The intersection is regarded as 
the secondary vertex. Figure 1 shows the scatter plot 
of the 7r + 7r _ invariant mass versus the decay length 
in the transverse plane (L xy ) for candidate events. 
The cluster of events with mass consistent with the 
nominal Kg mass indicates a clear Kg signal. Re- 
quirements \m^-„- - 0.497| < 0.018 GeV/c 2 and 
L xy > 0.01 m are used to remove background from 
non-Kg events. 
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Fig. 1. The distribution of m^^ versus L xy for oppositely 
charged track pairs in 7r + 7r~ K^-k^ candidate events. 



Events are kinematically fitted with four constraints 
(4C) to the hypothesis ip(2S) — * ir^n^ K ± tt t . Here 
^1^2 are associated with the K s decay, as determined 
above. For the remaining two tracks, the identifica- 
tion as K or it is done in the following way: (i) If the 
momentum of one track is larger than 1.34 GeV/c, 
that track is assigned to be a kaon and the other track 
a pion. (ii) If the momenta are both less than 1.34 
GeV/c, the fit is applied to the two possible combina- 
tions, and the one with the smaller x 2 is chosen. The 
confidence level of the selected 4C fit is required to be 
larger than 0.01. 

In addition, the combined chisquare (Xcom) f° r 
the assignment ip(2S) — > ir + Tr^K ± Tr^ is required to 
be smaller than those for the alternative hypotheses 
ijj(2S) -> 7T+7T-K+K- and ijj(2S) -> 2(tt + tt-). Here, 
the combined chisquare, Xcom> i s defined as the sum 
of the x 2 values of the kinematic fit (xline) anc ^ those 
from the particle identification assignments of the 

four tracks (xpid) <E1 : xL» = E l XjW*) + xL- 
After the above selection, the Dalitz plot for 
KgK^n 11 candidate events, shown in Fig. 2(a), 
is obtained. Monte Carlo simulated ip(2S) — > 
K*(892)K + c.c. events, shown in grey-scale, lie along 
a horizontal band for K*(892)°K + c.c. events and 
a vertical band for K*(892) + K~ + c.c. events. Fig- 
ure 2(b) shows the K ±r K T invariant mass after an ad- 
ditional requirement m K o^± > 1.0 GeV/c 2 to remove 
K*(892)+K- + c.c. events, and Fig. 2(c) shows the 
Kgir^ invariant mass after an additional requirement 
m K ±^ > 1.0 GeV/c 2 to remove K*(892)°K° + c.c. 
events. 

Contamination from background channels, which 
pass the selection criteria for the K^K^tt 11 events, 
mainly come from ip(2S) — > "fXcj, XcJ -^KgK + ir~ + 
c.c, which are assumed to take place via the inter- 
mediate state K*(892)K + c.c. Using the branching 
fractions taken from the PDG (Q), the contamination 
from these channels is estimated to be less than 0.6 
events for both the charged and neutral modes. The 
events along the third side of the Dalitz plot (diagonal 
side) in Fig. 2(a) may be due to the process ip(2S) — > 
p(2150) ± 7T = F j p(2150) ± -> K^K*. These contamina- 
tions are not corrected for, but are included in the 
systematic error, they and contribute an additional 
5.6% systematic error for both the charged and neu- 
tral modes. The contributions from other backgrounds 
are negligible. 

The invariant mass spectra for X ± 7r =F and Kg7r T 
are fitted using the if* (892) signal shape determined 
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Fig. 2. (a) Dalitz plot, (b) K^n^ invariant mass, and (c) 
Kgir^ invariant mass for KgK^n^ candidate events after the 
selection described in the text. The circles are data, and the 
shaded regions are simulated ip{2S) — * K* (892)K + c.c. events 
(the horizontal band for K* (892)° K + c.c. and the vertical 
band for K* (892) + K~ +c.c). The curves in (b) and (c) show 
the fit described in the text. 



with MC simulation plus a second order polynomial 
background and an additional Breit-Wigner func- 
tion for the if} (1430) (described below), as shown in 
Figs. 2(b) and (c); 65.6 ± 9.0 K*(892)°K° + c.c. and 
9.6 ± 4.2 K*(892) + K~ + c.c. events are observed. 
Their detection efficiencies are (9.68 ± 0.07)% and 
(7.25 ± 0.07)%, and their statistical significances are 
ller and 3.5a, respectively. 

In addition, 10.5 ± 5.1 events and 11.2 ± 5.3 events 
near 1430 MeV are found in the invariant mass spec- 
tra of K ± tt t and KgTT ± , respectively, by fitting with 
Breit-Wigner functions with the means fixed at 1.43 
GeV/c 2 . Their fitted widths are roughly 46 MeV/c 2 
and 100 MeV/c 2 , and corresponding statistical signif- 
icances are 3.4cr and 3.1<r, respectively. These events 
might be associated with the X *(1430), if (1410) or 
^(1430), but the limited statistics does not allow a 



determination of the spin J (=0, 1, or 2). Their de- 
tection efficiencies are (9.2 ± 1.0)% and (7.7 ± 0.9)%, 
where the errors reflect the unknown spin. 



3. Systematic errors and contributions from 
continuum 

The branching fraction is calculated from 

-,obs 



B 



fc 



where n obs is the number of observed K* (892)K + c.c. 
events, e is the detection efficiency obtained from the 
MC simulation, f c — (96. 3 ±3. 3)% is an efficiency cor- 



rection factor from the Ks reconstruction Ijllfl . N^^s) 
is the total number of i/j(2S) events, and Bi nt = 1/3 is 
taken as the branching fraction for K* (892)K+c.c. — > 
KsK^-k^. The K s — ► 7T + 7r branching ratio was in- 
cluded in the Monte Carlo simulation. 

Many sources of systematic error are considered. 
Those associated with the efficiency are determined by 
comparing J ftp and ip(2S) data with Monte Carlo sim- 
ulations for very clean decay channels, such as J/tp — > 
pn, K* {892) K + c.c, and ip(2S) -> tt+tt" J/ip, which 
allow the determination of systematic errors associ- 
ated with the MDC tracking efficiency, kinematic fit- 
ting, particle identification, photon selection efficiency 
and other experimental effects l|l3|) . The uncertain- 
ties of the background shapes and the total number 
of ip(2S) events are also sources of systematic errors. 
The total systematic errors for charged and neutral 
K*(892)K + c.c. mode are 14.0% and 12.6%, respec- 
tively. Table 1 summarizes the systematic errors. 

Contributions from the continuum e + e~ — ► 7* — > 
hadrons are estimated using a sample taken 

at = 3.65 GeV of 6.42 ± 0.24 pb" 1 Q, about 
one-third of the integrated luminosity at the ip{2S). 
In K*{892)°K° + c.c, 2.5+f;| events are observed, 
as shown in Fig. 3, while no events are observed in 
the K* (892) + K~ + c.c. channel, which corresponds to 
O.Otj^ events at the 68.3% confidence level jl7^ - 



4. Results 

Table 2 summarizes the observed numbers of 
events, detection efficiencies, and corresponding cross 
sections for K*{892)K + c.c. channels at y/s = 3.65 



GeV and wim^S)-. 
ficiencies at y/s ■■ 



respectively. Here, detection ef- 
= 3.65 GeV include the effect of 
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Table 1 

Summary of systematic errors (%). 







is* { Qa r >\0 IS® | ~ „ 
A ^oyzj J\ -\- c.c. 


Tracking 


8.0 


8.0 


kine. fit 


2.9 


2.9 


Bkgd shape 


5.0 


5.1 


Bkgd comtam. 


8.4 


5.7 


reconstion 


3.4 


3.4 


MC statistics 


1.1 


1.0 


N ip(2S) 


4.0 


sum 


14.0 


12.6 



initial state radiation 
tions of e + e 
and e + e~ — 



In addition, 
K*j(U30)"K° + c.c. 



K*j(U30)+K- 



c.c. 



the cross sec- 
ifgif ± 7T = F at 

3.686 GeV are determined to be 6.0 ± 2.9 ± 1.0 pb 
and 7.7±3.6±1.3pb, respectively, where the first 
errors are statistical and the second systematic, in- 
cluding the uncertainty from the unknown spin. The 
corresponding upper limits are 11.4 pb and 13.7 pb, 
respectively, at the 90% confidence level. 

Table 3 lists the branching fractions for the 



10 



C\J 
03 6 

C5 

"5" 4 

CM 

E 2 


CM 

to 

2 



(a) - 



10 



m 2 (K 7i) (GeV/c 2 ) 2 



5 



> 



- 1 — i — i — I — i — i — i — i — I — i — i — i — r 



(b) 



L 
0.6 



I I . .llll.Hl 



1.1 1.6 2.1 

m(Kjt) (GeV/c 2 ) 



2.6 



- 1 — I — I — | — I — I — I — I — | — I — I — I — I — | — I - 

(c) 



_l I I I I L 



J I I I I L 



0.6 



1 .1 



1.6 



2.1 



2.6 



Fie 



m(K s 7t) (GeV/cT) 

3. (a) Dalitz plot, (b) K^ir^ 1 invariant mass, and (c) 



K'gTT^ invariant mass for K^K^tt^ candidate events in the 
\/s = 3.65 GeV data sample. The circles are data, and the 
shaded regions are simulated ip(2S) — > K* (892)K + c.c. events 
(the horizontal band for K* (892)° K + c.c. events and the 
vertical band for K* (892) + K~ + c.c. events). 



Table 2 

Observed cross sections for e + e — — > K* (892)K + c.c. at y/s ■■ 
3.65 GeV and 3.686 GeV. 



Channels 


v 7 ^ 
(GeV) 


jqobs 


e 
(%) 


(7 

(Pb) 


K*(892)+K- + c.c. 
K*(892)°K° +c.c. 


3.65 


0+ 1 ' 3 

o c+2.6 
z -°-1.8 


5.7 
7.7 


< 21 (90% C.L.) 

16±\l ± 2 

< 42 (90% C.L.) 


A"* (892)+ A"- + c.c. 
K*(892)°K° +c.c. 


3.686 


9.6 ±4.2 
65.6 ±9.0 


7.3 
9.7 


20.9 ± 9.1 ± 2.9 
107 ± 15 ± 13 



Table 3 

Branching fractions measured for ip(2S) — > K*(892)K + c.c. 
The corresponding J/tp branching fractions (2) and the ratios 



h B(J/V) 


also 


given. 






Channels 


BW>(2S)) 
(xl0~ 5 ) 


B(J/V0 
(xlO- 4 ) 


Qh 
(%) 


K*(892) + K~ H 
K"*(892)°A° + 


- c.c. 
c.c. 


2-9±i;| ± 0.4 
13.3+Jg ± 1.7 


50 ± 4 
42 ±4 


59+ ' 27 
3.2 ±0.8 



^(25) -> AT* (892).^ + c.c. decay modes, where the 
contributions of the continuum is subtracted in- 
coherently with normalized integrated luminosity 
without considering its interference with the reso- 
nant amplitude. The table also lists the branching 
fractions of J/ip decays (0) as well as the ratios of 
the tp(2S) to J ftp branching fractions. The ratio 



B(ip(2S)->K* (892)" K +c.c.) _ 



= 4.6 



-2.9 



B(v(2S)^if* (892)+ k- +c.c.) — "*-"-2.2 snows a large 
isospin violation between the charged and neutral 

modes ofip(2S) -> K*(892)K + c.c. decays. Since the 
amplitudes for ip(2S) — > if*(892)X + c.c. decays up 
to first order in the SU(3) breaking consists of two 
parts: the strong de cay amplitude and the electro- 
magnetic amplitude l|18|) . a possible interpretation for 
this large isospin violation is that the electromagnetic 
decay amplitude plays an important role in "0(25) — > 
K*(892)K + c.c. decays, while in J/tp decays the 
strong decay amplitude dominates. The results listed 
in Table 3 show that, ip(2S) -> K* (892)1? + c.c. 
branching fractions are strongly suppressed with re- 
spect to the pQCD expectation. The charged branch- 
ing fraction is suppressed more than the neutral one 
and is consistent with the upper limit measured by the 
Markll experiment (< 5.4 x 10" 5 , at 90% C.L.) (Q). 
Our results marginally accommodate the predictions 
by Chaichian et al. and by Feldmann et al. (5), two 
predictions for the charged mode branching fractions 
are 4.5 x 10~ 5 and 1.2 x 10~ 5 respectively, while being 
larger than their two predictions for the neutral mode 
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(7.6 x 1CT 5 and 5.1 x 1CT 5 , respectively). Ma (g) pre- 
dicted Qk-k to be (3.6 ± 0.6)%, which is consistent 
with our measurement for the neutral mode, but not 
for the charged mode. 

Based on the observed cross sections in Table 2, the 
branching fractions of ip(2S) -> K*(892) + R- + c.c. 
and K*(892)°~K° + c.c. may be calculated by the 
model proposed in Ref. (|l5T l. where all the contribu- 
tions from the continuum one-photon annihilation 
amplitude, the electromagnetic amplitude and the 
three-gluon amplitude of the il>(2S) decay, and their 
interferences are taken into account. By fitting these 
observed cross sections the phase between the elec- 
tromagnetic amplitude and three-gluon amplitude of 
the V , (25 l ) decay is constrained in the range from 95° 
to 304°, disfavors the positive solution of the orthog- 
onal phase ±90° determined from J/ip decays (|6j). 
The branching fractions in this case are: B(ip(2S) — > 
K*(892)+K- + c.c.) = (3.1+iJ) x 1CT 5 and 
B(<ip(2S) -» K* (892)° K° + c.c.) = (13.7±^;g) x 10" 5 , 
where the large errors are due to the large phase 
uncertainty. 

In conclusion, we present the branching fractions 
for tp(2S) -> K*(892)°K° + c.c. and K*(892)+K- + 
c.c. They are suppressed with respect to the pQCD 
expectation, and a large isospin violation between the 
charged and neutral mode is observed. These results 
are compatible with those recently reported by CLEO 
Collaboration for the same channels in ref. l|19l) . 
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